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We investigate, in the case of the 2p — Is transition in atomic Hydrogen, the behaviour of the 
spontaneously emitted electromagnetic field in spacetime. We focus on Glauber’s wave function for 
the emitted photon, a quantity which we find is nonzero outside the lightcone at all times after 
the start of the emission. We identify the uncertainty on the position of the decaying electron 
as a source of departure from causality in the naive sense of the term. We carry out a detailed 
study of the emitted electric field in the mid- and far-held regions, through analytical and numerical 
computations as well as asymptotic arguments. 


I. INTRODUCTION 

Fermi was the first to study the electromagnetic field 
emitted during an atomic transition, which has since 
then been a recurring theme of investigation in atomic 
physics and quantum electrodyamics [IHi]. Fermi [T] 
found that, if the survival probability of the electron 
in the excited state is assumed to decay exponentially, 
according to the usual Wigner-Weisskopf aproximation 
[H] , then the emitted electromagnetic field will propagate 
causally, in other words, vanish outside the lightcone 
centred around t — 0 (the instant at which the emis¬ 
sion starts) and the position x = 0 of the Hydrogen 
nucleus (proton). It was then noticed by Shirokov 
[2] that Fermi’s result made use of an unmentioned 
approximation, which consisted in extending the range 
of integration over electromagnetic frequencies from 
the positive real semi-axis to the whole real axis, 
thereby including nonphysical electromagnetic negative 
frequency modes in the treatment. Hegerfeldt later [?] 
formalised and generalised Shirokov’s remarks, linking 
the absence of negative electromagnetic frequencies with 
the noncausal field propagation via the Paley-Wiener 
theorem [7] on holomorphic Fourier transforms. 

Note that in [IHl] the authors consider not one atom 
in free space but two, the hrst one being initially in its 
excited state while the second one starts in its ground 
state. Rather than computing the emitted field, they 
focus on the probability of excitation of the second atom 
as a function of time (this is known in the literature as the 
“Fermi problem”). Our present study features a single 
atom, initially in its excited state, and purports to obtain 
the spacetime dependence of the spontaneously emitted 
field, expressed by Glauber’s photon wave function. 

Our main new result is the following: even if we 
artificially integrate over negative frequencies, thereby 
bypassing Hegerfeldt’s objections, causality is still 
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violated, which establishes that Hegerfeldt’s mechanism 
is not the single source of noncausality in the present 
problem. This result follows from a rigorous computa¬ 
tion of the single photon wave function, and by the use 
of exact expressions for the coupling coefficients between 
electronic Is and 2p states of the Hydrogen atom and 
electromagnetic modes mm- 

In sect. [IT] we review the tools for the description of 
spontaneous emission. Sect. [m] sets the stage for the 
computation of the wave function of the emitted photon, 
which is the object which we use in order to assess caus¬ 
ality. In sect. |IV| we review the usual treatment Dun] 
where causality is derived by the way of multiple approx¬ 
imations (dipole approximation, extension of the electro¬ 
magnetic spectrum to negative frequencies, usual E • x 
coupling instead of the minimal A • p coupling). In 
sects. |V| and |VI| we refine the treatment by progressively 
waiving various approximations, and we find that the res¬ 
ult is no longer causal. We identify the different sources 
of noncausality, and study the spacetime dependence of 
the photon wave function in the mid- and far-held regions 
in detail. Sect. ivn is a discussion of our results. 

II. THE DECAY OF A TWO-LEVEL ATOM 

Let us consider a two-level atom (ground state | g), 
excited state | e)) interacting with the electromagnetic 
held. The atom sits in free space. The Hamiltonian H = 
Ha + Hr+Hj is a sum of three terms: the atom Hamilto¬ 
nian Ha, the electromagnetic held Hamiltonian Hr, and 
the interaction Hamiltonian Hj. In the Schrodinger pic¬ 
ture these read m 

TLa = |g)(g|+^e |e)(e| (la) 

Hr='^ j dfc fic 11 k| I (k,t = 0) 0 (a) (k,t = 0) 

Hi = —A (x,0) • p 
me 


(lb) 

(l c) 




2 


where Wg/e is the angular frequency of the ground/excited 
atomic level, A labels the polarisation of the electromag¬ 
netic field, rrie is the electron mass. Also, x is the electron 
position operator and p is the electron linear momentum 
operator. The field operator A (x, t) is the vector po¬ 
tential (here we work in the Coulomb gauge), which is 
expanded over plane waves as 


is the usual volume element on the lightcone [H HI 
(where 0 stands for the Heaviside distribution), which is 
an invariant under Poincare transformations. The polar¬ 
isation vectors (k) are any two mutually ortho¬ 

gonal unit vectors taken in the plane orthogonal to the 
wave vector k. Finally, we give the commutation relation 
between the photon ladder operators, which reads 


A (x, t) = X! / 

+a(A) (k, i) e(A) (k) . 

( 2 ) 

Here 

(27r) 


«(>.) (k), a|'^) (q) =2||k||(27r)^(5(k-q)5^A. (4) 


We consider spontaneous emission in a vacuum: at t = 0, 
the electron is in its excited state, while no photons are 
present in the field. For such an initial condition, the 
rotating wave approximation holds |14j and the state of 
the system at time t > 0 reads 


IV'(t)) =Ce(t)e--"‘ |e,0)+ ^ [ d'A:cg,A(k,f)e-'(-«+^ll'^ll)‘ |g,lA,k) (5) 

A=±-^ 


where | e, 0) means that the atom is in its excited state 
and the field contains no photons and | g, 1a, k) means 
that the atom is in its ground state and the field contains 
a photon of wave vector k and polarisation A. 


and I g, lA.k) ■ These are well-known for the 2p — Is hydro¬ 
gen transition. Writing, for this transition, | g) =| Is) and 
I e) =12pm2), with m 2 the magnetic quantum number of 
the 2p sublevel considered, one has [9] 


Let us turn to the matrix elements of the interaction 
Hamiltonian in the Hilbert (sub)space spanned by |e, 0) 


(Is, lA,k |lL/| 2 pm 2 , 0 ) 


h he 2^ 




(k) • I. 


CqC TTIq Q-q 


34 


1 + (§ao ||k| 


( 6 ) 


where we introduced the Bohr radius oq. The ^^2 are 
given by 


^0 = Gz, 


l±i = T 


Gx A i^y 


(7a) 

(7b) 


Since we are interested in spontaneous emission, we set 
Ce (t = 0) = 1 and VA G {l,2}Vk G R^Cg,A (k,t = 0) = 0. 
We want to compute probability amplitudes of emission, 
namely 


In order to derive §, one must remember the expressions 
for the electronic wave functions of the Is and 2p m 2 sub- 
levels: 


Cg.A (k,t) = (g, lA.k|t/(t)|e,0) (9) 


V'ls (x) 


'02pm2 (x) 




x-l 


m 2 ■ 


(8a) 




(8b) where (7(t) = ex, ^ 

for the system. From Q and ([^, we get 
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Ce it)=-J 2 


dk 


A=± ' 


(2^r2||k|| 


Gl (k) Cg,A (k, t) e-i(d|k||-<^e+c^g)t^ 


Cg.A (k, t) = --Gx (k) Ce (i) 


where Gx (k) = (Is, lA,k I-&/I 2p m 2 , 0). 


quantum optics. 


(10a) 

(10b) 


Consider a pure, single-photon state of the electromag- 

III. THE PHOTON WAVE FUNCTION netic field 


The single-photon wave function, hereafter referred to 
as “the photon wave function”, is a very useful object 
for the description of one-photon Fock states of the elec¬ 
tromagnetic field either in momentum space m or, as 
developed more recently [51IIMIS1, in direct space (and 
time). In the rest of the paper, we will use Glauber’s 
photon wave function to investigate the spontaneous 
emission of light during the atomic transition at hand, 
and focus on causality. This wave function was first 
introduced by Glauber and Titulaer |19j in order to 
characterize correlations of the electromagnetic field in 


|l:/(i)) = Xl / (k,0a(A) (k) |0). (11) 

The photon wave function can be obtained through 
Glauber’s extraction rule which, when states and operat¬ 
ors are defined in the Schrodinger picture, reads [IMI^ 

tAi(x,t) = (0|Ei(x,0) |l,/(t)) (12) 

where Ej^ represents the transverse part of the electric 
field operator defined through 


A=± ' 


Ex (x,t) = iW—y] / dfc||k|| a(A) (k, t) e(A) (k) (k, t) (k) e 


....—ikx 


(13) 


In our problem, the state of the electromagnetic field is not pure, but rather, as seen from ([^, entangled with 
that of the atom. Hence, projecting onto the single photon sector and applying Glauber’s extraction rule (12), the 
single-photon wave function reads 


i>± (x, t) = (g (t), 0 1 Ex (x, 0) I i/j (t)) 

= / d/c (01 Ex (x, 0) Cg,A (k, t) | 1a, k) 

A=± 

= n/r E / dfc||k||e‘(‘^--ll‘^ll‘)cg,A(k,t)e(A)(k) (14) 


A. Formal computation: preliminary steps 


At this point, it comes in handy to notice that (10b) can be formally integrated, yielding 

Cg,A (k, dt' Ga (k) Ce {t') eHdIkI I 

so that the single-photon wave function reads, in the most general case of our problem, 

^x (X,t) = J dfc||k||e’d^'"-dlk|l‘)e(^) (k) Ga (k) dt'Ce (t') 


A=±- 

. 25 He 
34 eomeOo 


5] / dfc||k||e‘d^— 


<Ai (k) ■ 


A=±' 


m 


dt'Ce (t') 


(15) 


(16) 


I 

where we introduced luq = uig — Wg and fcx = 3/ (2ao). The unit polarisation vectors obey the closure relation 


EGw) (k)e{,)(k)=5*^ 


A=± 


k2 


(17) 
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so that the wave function now is 


l/)_L (X, t) 


.22 he 
34 eoTOeao 


dkllkl 


gi(k-x-c||k||t) 







(18) 


Choosing a coordinate system for which x points along the third axis e^, we can compute the angular integrals: 

r27r 


pTT pZTT 

F(fc, ||x||)= / d0sin0 / di^ 
Jo Jo 


^m2 1 


k \ k 

■ l|k||j ||k|| 


d^ll 


Ifc X COS t 


= d6»sin6»e*'^'ll’‘ll™"® / ■ 

Jo Jo 


= 27r/ d6»sin6»e**^ll’‘ll™"® 


.27r 


■ p{x) ■ 

Sm2 


sin 9 cos ip 

/ dv5< 


to 

- 

sin 9 sin p 

/o 


AA 


cos 9 


1 




(^m] sin dcos(p + sin 6* sin </? + cos 9 


r^O 


(1 - I sin^ 9) 


CL"i(l-isin2 0) 
sin^ 9 


= 27rl(fc, ||x||). 
The integrals over 9 give 


&’y'> {k, 


X = 1- 


A^^y) 

: S 1 TI 2 


fc llxl 




-ifc||x|| _ gifcih 


(fc,||x||) = -2i 


A^) 

• Sm2 


fc ||x|| V 

_ )_ ( e”* 

k ||x|| V 


„ifc||x||\ ^ /"p-ifcllxll _pifc||x||\ 

/ (fc||xin"v ) 


-ife||x|| _|_ gifc||x| 


(fc||x||)' 




-ifc||x|| _ gife||x| 


(19) 


( 20 ) 


As we can see, the photon wave function contains contributions proportional to 1/ ||x||, 1/ ||x|| and 1/ ||x|| , which 
are respectively known as the far-held, mid-held and near-held contributions. We then have 


i>± (x,t) = -i 


23 he f+°° dk 


-.—ickt 


34 eomeOo Jo (27r)" 


1+ £7 


yj) 


^I(fc,||x||) /‘dt'ce(Oe'(“'=-^“^*'. (21) 

Jo 


In most of what follows we will use the Wigner-Weisskopf 
approximation of exponential decay. We shall then have 


Ce(t) =e-‘‘"^®*e-3r* (22) 


where wls is the partial Lamb shift [20] of the excited 
2p level due to the Is level and T is the decay rate. This 


yields 


pt pt 

/ dt'Ce (t') = / dt'e*(^'=-(‘^'>+‘^^sii*'e-3rt' 

Jo Jo 


gi(c/c-(cJo+^Ls))i' 
i(cfc- (wo -hWLs)) - 
X _ gi(c/c-(wo+WLs))t g-jrt 
ck- (wo +wls) + Ar 


We introduce the space-saving notation. 


(23) 


= Wo -f wls ~ 2^ I 


(24) 
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We can then write the contributions to the far-held, mid- held, and near-held to this photon wave function; 


J.,y) .2i he dfc e- 

-Ti—Ti- / - -k- 


-L(far) 34 eometto ||x|| Jq ( 27 r)^ 


V' 


_ 2i he /■+“ dfc 


' + (>&)■ 

^ —iflnt 


1 2 






-l-(mid) 34 epTOeOo ||x||^ Jo (27r) 


, ,212 (® ' 


-ifc||x|| _(_ gife||x| 


= -1 


.25 he /■+°° dfc 1 e-^^o* 


_L(near) 34 ||x ||3 k 


1 + 


(*)■ 


-i'=l|x|| _ gifc||x| 


V'Sfar) = 0i 


V' 


( 2 ) __ 

±(mid) ^4 


2i fie e 


(z) g-l -00 


dfc 


eouieO-o ||x||^ 7 o ( 27 r)^ '' 


—iQo^ 


n 2 


■0 


(z) _.2i fie Cml /■+“ dfc 1 e 




^g-ife||x|| _|_ gifc||> 


^ _ g— i(cfe —f 2 o)t 
ck — VLq ’ 

2^ _ g-i(c/c-no)t 

cfc — rio ’ 

2 _ g-i(cfe-no)t 


cfc — flo 


-L(near) 3^ eoWeOo | |x| io (27r)^ ^ 


1 + 


(0 


1 2 




-ife||x|| _ gifc||x| 


2 _ Q-i(ck-flo)t 

ck — flo ’ 

2 _ g-i(cfe-no)t 
cfc — rio 


(25a) 

(25b) 

(25c) 

(25d) 

(25e) 

(25f) 


B. General method 


We write 


^- 1-00 

/ dfc-- 
Jo 


'2—n 


1+ £7 


(fi)' 


^+00 

/ dfc/„(fc)e±*'=ll’ 

Jo 


n 2 ' 


±ifc||x| 


2 _ Q-i{ck-ujo)t 

ck — ujQ 


It is a general result of distribution theory m that 


(l|x|| ,0 = 


(5(-)- -Vp- ) 

TT 


(T| 


= lfn{T\\x\\,t)- ^ lim 

2 27r £->-0+ 


^+00' 


da 


T||x|| -u 


fn (cr, t) 


(26) 


(27) 


with the Fourier transform 


IV. CAUSALITY IN THE STANDARD 
TREATMENT 


^ + 00 


(28) 


/„(||x||,t)= / dfc/„(fc,t)e fu f 11 • r u u m rmi 

J_^ Here we use the following standard [IJ [lOj approx¬ 

imations used to established the causality of the wave 
function of the emitted photon: the Wigner-Weisskopf 
In what follows we will endeavour to compute the exponential decay, the usual E • x coupling between the 

integrals, hrst by making use of several aproximations in held and the atom instead of the minimal A • p coupling, 

order to obtain a causal result, and thereafter progress- the dipole approximation and the approximation which 


ively waiving these approximations. 


consists, as we will see, in extending the range of elec- 
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Figure 1: Jordan loops in the complex fc-plane used to 
compute the integrals (30). The (isolated) simple pole 
rio/c of the integrands is represented by a red circled 
cross. 


tromagnetic frequencies to the negative real semi-axis. 


We switch (only in the present section) from the min¬ 
imal A • p to the usual E • x coupling. In the dipole ap¬ 
proximation, which, for the A-p coupling, consists [22] in 

\2 


1 -h {k/k^Y 


denominator in (21), 


forgetting about the 

this substitution results [H] [M] in the multiplication of 
the interaction matrix element by c/c/wq. Plugging this 
in (21), we have 


■0_L (x,t) 


7 

23 hec 
eoTUettoUJo 



(2^)2 


l|x||) 


g—i(cfc—Ho)* 

ck — fig 


(29) 


The usual trick [Diiiin] is then to extend the range of 
integration from the positive real semi-axis to the whole 
real axis, and we find ourselves computing integrals of 
the type 


hSL) (II^i 


,t) = 


^+oo 


dk 


2^ _ g-i(cfe-no)t 

ck — fig 


(30) 

with n G {1,2,3}, and where the label -!-(—) is assigned 
to outgoing (ingoing) radial waves. Here the std sub¬ 
script stands for “standard” as we follow the lines of the 
standard treatment Dun] of the problem. For this we 
use Cauchy’s residue theorem (see Fig. [^. Taking into 
account the fact that t and ||x|| are positive quantities, 
we fi nd ( for more details, see the similar treatment of 
sect. 


VA) 


<Id)(IWI,t)=2|0(ct-||x||)(^^y "eWciixii 

<sld)(IWhi)=0. 

(31) 

This means that contributions from ingoing waves are 
zero, as found for instance in m- As made clear by the 
Heaviside step, this result is explicitly causal. We finally 
have 


Y’± (x,t) 


23 he 

3‘^7r egTOeOgWo 


6»(ct- ||x||) 





(^IWI)^ (^l|x||)^ 

I i 1 


(ai||x||)^ (eailxll) 




(32) 


V. MINIMAL COUPLING IN THE DIPOLE previous sect. |TV) This yields 

APPROXIMATION 


A. Intrinsic noncausality in the near-fleld 


We now return to what we regard as the more correct 
coupling: the minimal A • p coupling. We still work in 
the dipole approximation and in the Wigner-Weisskopf 
approximation. We shall first extend the range of integ¬ 
ration to the negative real semi-axis as was done in the 


(x,i) = 


22 he 
3^ egmeOg 



dfc /;2g-inot 
(27r)^ 


l|x||) 


g—i(cfe—no)t 

ck — Hg 


(33) 
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Figure 2: Jordan loops in the complex fc-plane used to 
compute the integrals (34). The (isolated) simple poles 
rio/c and —\a of the integrands are represented by a red 
and a cyan circled cross, respectively. 


We find ourselves computing the following Fourier trans¬ 
forms: 


/n(dip) (T 11^11)0 — 


^+oo 


d/cA:2-"e=^ife|l> 


^ _ Q—i{ck—Qo)i 


' —oo 
/*+oo 


ck — rio 

dA/„(dip) (A:)e±‘'=IWI. 


(34) 

Here the dip subscript stands for “dipole” as we use the 
dipole-approximated minimal A • p coupling. We use 
Cauchy’s residue theorem (see Fig. [^. Taking into ac¬ 
count the fact that t and ||x|| are positive quantities, we 
find, for n G {1,2} 


/n(dip) (^; 0 

2i7r 



2 —n 


[—9 (m) -I- 6* (u -I- ct)] 


(35a) 


whence 

/n(dip) (||x|| ,t) = 0 (35b) 

and 


/n(dip) (- I |x| I , 1) = 2i7r 




+ ct ), 


from which we conclude 


/n(dip) (- l|x|| ,t) T /n(dip) (l|x|| ,t) = 2i7r 



{?l|x||0(_||x||+ct) 


(35c) 


(35d) 


r 


which is quite transparently causal (notice from (35b) 
that the contributions from incoming waves are identic¬ 
ally zero here, as they were in the standard treatment of 
sect. IV). Now, notice that for n = 3, the integrand in 
ras a (simple) pole at fc = 0 (see Fig. [^. Since we 


(34) 


integrate over the real axis, this could be a serious prob¬ 
lem, because it means that are divergent integrals. 
We are saved from dealing with such divergences by noti¬ 


cing that, according to (|19|), we will only be interested in 

|x|| ,t), which 


computing the difference — Jg 

is an integral over a function which, according to (34), 
has only an artificial singularity at fc = 0 [53]. Accord¬ 
ingly, we can shift this singularity away from the real 
axis 0 ——ia to compute the integral (see Fig. [^, be¬ 
fore taking the limit a —>■ 0 at the end. This amounts 
to performing the substitution -G {k + ia) " and 
yields 


/3(dip) (- l|x|| ,t) - /3(dip) (||x|| ,t) 


2i7r 


c 

Hg L 


0 {ct — 




if2ol|x|| 


- + (l - 


(36) 


Comparison with (35d) shows that the result for n = 3 
features not only an extra causal contribution, but also 


a completely noncausal term. This is a feature of the 
slow dependence of the A • p atom-field coupling on the 




























norm of the electromagnetic wave vector k. As far as 
we know, similar calculations [nisiEiiio] of the outgoing 
field have mostly been carried out in the Power-Zineau- 
Woolley picture m of quantum electrodynamics where 
the interaction Hamiltonian is of the usual E • x form. In 
this case no pole is present at A: = 0 (see Figs. and[^ 
and one retrieves a causal result, as was done in sect.lIVl 


B. Hegerfeldt theorem for the minimal coupling in 
the dipole approximation 

We now waive the approximation which consists in ex¬ 
tending the range of integration to the negative real semi¬ 


axis as was done in the previous sect. IV This yields 


(x,f) = - 


22 he 


^+oo 


34 


eomeflo Jo 


dk 


iOot 


1 _ p-i(cfe-no)t 

I(fc.lNI)' o (37) 


ck — rio 

and we find ourselves computing integrals of the type 


<d\p) (11^1 


r+o° 1 p—i(cfc—O q)* 

\,t)= 

Jo 


I 

J —c 


ck — Ho 

dA:0(fc)/„(dip) (fc)e±*^-IWI 


with n S {1,2,3}. With (27) in mind we can use 


the Fourier transform /n(dip) of /„(dip) computed in 
sect. ra Since we did not extend the range of in¬ 
tegration to the whole real axis, we will not retrieve 
a causal propagation, as first noted by Shirokov [2]. 
Note that this is true regardless of the choice for the 
coupling: Fermi’s proof of causality used the usual E • x 
coupling, and included, as an approximation necessary to 
causality, the extension of the range of integration to the 
whole real axis. This is an illustration of the Hegerfeldt 
theorem [3] , which states that non-causalities will always 
arise for Hamiltonians bounded by below. The relevant 
Hamiltonian here for the Hegerfeldt theorem is the 
Hamiltonian Hr ([^ of the free field, which has K+ as its 
spectrum and is hence bounded. The non-causality seen 
in (36) is of a different kind—it is not a manifestation 


of the Hegerfeldt theorem—and is not known in the 
literature. See sect. m for further discussion. Notice 
that the n = 3 case at hand corresponds, as can be seen 
from (19), to the near-held part of the emitted photon 


wave function. In the present manuscript we will rather 
focus on the mid- and far-held contributions to the 
electric held, and shall return to the delicate question of 
the near-held in an upcoming manuscript. 

We now focus on the Hegerfeldt noncausality. To in¬ 
vestigate this particular point, we compute the convo- 


(38) (ne{l,2}) 


lutions of (|35d[) and (|36|) with the principal value as 
prescribed 


Dy (27). In the mid- and far-held we get 


C'n(dip) (I |x| I , t) — 


27r 


Vp - * /n(dip) {-,0 


= e 


Hq 

c 


2—n 


p. ( .^0 II I 

—El I 1— X 
c 


Ei ( i— (||x|| -I- ct) 


(39) 


and in the near-held 


Csidip) (l|x|| ,t) = = - 


= e 


ill 


27r 


-vp 


vp - */3(dip) i-.t) 


lA 

Hn 


-Ei(i^||x| 


-Ei ( i^ (||x| 


ct) 


l• + 00 

^ _ gifiot 

r+oo 

du 

y-iixii 

V J 

-||x||- 

-ct u 


(40) 


SO that 


^3(dip) ( 11^11 ^3(dip) 


= 


C 

iln 


—e 


_log: , 

Ho °Hi-iixii+cfi; ■ 


-Ei ( -i^ I |x| 1^ -I- Ei ^i^ (- I |x| I -I- ct) 
-Ei ( i^ ||x||^ + Ei (i^ (||x|| -hct) 


ct 


(41) 


Here Ei stands for the exponential integral 

e“ 


Ei (x) = — 


du¬ 


ll 


(42) 


Note that to deduce /„(dip), as given by (35b), from 
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(27) and (^34|, we made unwarranted use of Jordan’s 


lemma: the integral of /„(dip) over the semicircle from 
7top (see Fig. is not zero for n = 1 in the limit of large 
semicircle radius. On the contrary, this integral diverges 
in that limit [IS]. It is thus clear that the dipole ap¬ 
proximation forbids a clear, consistent treatment of the 
problem at hand (unless a cutoff is introduced around 
the frequency wx = (3c) / (2ao) as explained in [HI [^ . 
However if a cutoff is introduced we cannot avoid 
Hegerfeldt-type non-causalities). In the next sect. VI 
we do away with this approximation. Nevertheless, we 
will see that very similar terms to those obtained here 
in the framework of the dipole approximation arise. 
Since the expressions obtained are much more involved 
in the next sect. lYll we study the less complicated 
results of the present section in some detail, and this 
knowledge will come in handy for later. Namely, 
taking the fcx —> +oo limit of the results of sect. |VI| 
yields the results in the dipole approximation. This 
confirms that the exact coupling provides the correct 
regularisation for the dipole approximation (see mm)- 

We plot in Fig. |^the square moduli of the causal (or 
pole) (35d I and noncausal (or principal value) ( [^ con¬ 
tributions to the far- and mid-field parts of the emitted 
photon wave function. 

Two noticeable patterns emerge: 


Inside the lightcone, we notice that the contribu¬ 
tions to the far-held (n = 1) from the pole term 
(35d I and from the principal value term (39) are al¬ 
most indistinguishable, at least when ||x|| < c/^Iq. 
This feature can be explained by simple asymp¬ 
totic [24] arguments: under the condition that we 


pick a spacetime point reasonably “deep” within 
the lightcone, we may make the approximation that 
I |x| I / {ct) —?► 0 so that from (39) and the asymptotic 
series 


p“ h\ 

Ei(u) - — (43) 

n—>-+oo u " 




we may write, with the extra help of the Taylor 
series for Ei (u) around m = 0, 

Cqdip) (-||x|| ,t)-C'i(dip) (||x|| ,t) ^ - iTT^, (44) 

in agreement with the asymptotic behaviour of the 
causal part (35d) in the same (Hq/c) ||x|| <C 1 limit. 


This good agreement is not reached for the mid- 
field because the contributions from ingoing and 
outgoing waves are added instead of substracted, 
as seen in (19). 


Outside the lightcone (where the only contributions 
to the far-and mid-field, obviously, come from the 
principal value integrals), we notice on the graphs 
that the {n = 1) far-held decays not as l/||x|| 
as expected, but as l/||x|| . As for the (n = 2) 
mid-held, it decays not as 1/ ||x||^ as expected, but 
as l/||x||^. Again, this can be explained by the 
asymptotic behaviour |24j of the exponential in¬ 
tegral function: under the condition that we pick 
a spacetime point reasonably far away from the 
lightcone, we may make the approximation that 
||x||/(ct) —>■ -boo so that from (39) and (43) we 
may write the asymptotic series 


f^l(dip) ( l|x|| ,^) ^l(dip) (l|x|| , t) — 
C' 2 (dip) ( I |x| I , t) -|- C 2 (dip) (I |x| | , t) — 


This accounts for the spacewise-decay behaviour 
described just above. The conclusion reached is in¬ 
teresting: outside the lightcone, the far-held decays 
more strongly than usual with increasing distance 
and mimics the usual behaviour of the mid-held 
(l/||x||^), while the mid-held decays much more 
strongly than usual with increasing distance, so 
that its l/||x||^ decay is stronger than the usual 
1/ ||x|b decay of the near-held. 

Both these points are noteworthy features of the 

Hegerfeldt-noncausal terms, which come from the fact 


-I- 

^ 11 il -i (for large^lMI), 

l|x|| 

2 + (_2 + 2ifIof) , 

/n^N2,, ,,9 (^O.). 


(45a) 

(45b) 


that the integration is restricted to positive electromag¬ 
netic frequencies, as it should. The hrst point conhrms 
that the usual approximation consisting in extending the 
integration to the negative real semi-axis is fairly solid: 
inside the lightcone, we see (Eig. that the contribu¬ 
tion from the (principal value) noncausal terms (39) just 
about equals that of the (pole) causal terms (35d), which 


means that the result yielded by extending the range of 
integration would be sensible. As for what happens out¬ 
side the lightcone, we not only pointed the well-known 
fact that taking the absence of negative frequencies into 
account yields a nonzero result, but we noticed the in- 
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Figure 3: See rest of figure and caption below. 


teresting fact that the field decays more strongly with 
increasing distance than would be naively inferred from 
the usual behaviour of the mid- and far-field contribu¬ 
tions to the emitted electric field. Keep in mind that the 
noncausal contributions to the emitted field are small, as 
seen on Fig. (remember that the axes are logarithmic). 


VI. EXACT TREATMENT FOR THE MINIMAL 
COUPLING 

We now switch to a fully rigorous treatment of the 
problem. We use the minimal A • p form of the atom- 


field coupling. We no longer work in the dipole approx¬ 
imation, but shall instead use the exact interaction mat¬ 
rix element (§. We also no longer perform the Wigner- 
Weisskopf approximation but use perturbation theory at 
short times |26|- Finally, we do not extend the range 
of integration to the negative real semi-axis as was done 
in sects. IV and V A In a time-dependent perturbative 


treatment of the present problem, we approximate (21) 
to first order in time, which consists [55] in considering 
that Ce (t) = 1, so that the photon wave function is ap¬ 
proximated by 


ip± (x, t) 


22 he 
3"^ eometto 



dk 




e-'-«‘l(fc,||x||) 


2 ^ _ 

ck — loq 


(46) 
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X t 


F igure 3: Square moduli of the contributions to the far- (n = 1) and mid- (n = 2) fields from the poles as given by 
(35d) and from the principal value integrals as given by (39) as a function of distance from the nucleus, for fixed 
time in the left column, and as a function of tim e, for fixed distance in the right column. This is for the minimal 
A • p coupling in the dipole approximation (sect. VB|. Both axes are logarithmic on all figures. The boundary of 
the lightcone is signalled by a vertical line. The reader will notice the noteworthy fact that for fixed x = c/lxq^ the 
contributions from the mid-field and far-held are indistinguishable inside the lightcone. 


We hnd ourselves computing integrals of the type 

^+oo 




r 

Jo 


dk ■ 


1 + 


(*)■ 


n 2' 


p±ife||> 


I _ Q-i{ck-ujo)t 

ck — UJq 


(47) 


with n G {1,2,3}. The integrand is similar to that in 
(34), but features two extra poles at fc = ±ifcx, and its 
“Wigner-Weisskopf” pole at k = ujq/c sits on the real 
axis. Since the latter is only an artihcial singularity, we 
can shift it to the lower half plane as seen on Fig. As 


I 

was the case in the previous sect.|Vj for n = 3 the integ¬ 
rand has a singularity at fc = 0, which we are allowed to 
shift to A: = —ia since we will be interested in computing 


the difference (— H, 


t(-) 


, t), which is an integ¬ 


ral over a function which, as seen from (47), has only an 


artihcial singularity at k = 0. We can rewrite (47) as 


(l|x|| ,i) = 


0(k) (fc + ia)^-\ ^i,||,|| i-e-Hcfe-^oh 


1 + * 


(i)' 


ck — (wo — ic e) 


^+oo 


' —oo 
/•+00 


dke{k) /„(/c,t)e±‘'=ll’^ll 
dke{k) (5„ {k) - K 


(48) 


where it is implied that the limit a —)■ O'*", e —>■ 0+ should 
be taken outside the integral. Here the functions gn and 


hn read 


_ (fc -I- ia)^ " 1 



^ ck — {loq — ic e) 

[!+(*) J 



(49a) 
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Figure 4: Jordan loops in the complex fc-plane used to 
compute the Fourier transform (501. The (isolated) 
poles wq/c — ie and ±ifcx are represented by red circled 
crosses while the (isolated) simple pole —ia is 
represented by a cyan circled cross. 




2 —n 


^—i{ck—ojQ)t 


r, /1 

^ ck — (wo — ic e) 

[i+(it) j 



(49b) 


Once again we use (27) and hence need to compute the 
Fourier transform 


/n(l|x||,t) = 


^+oo 


dk fn (fc, t) e 


— \kx 


(50) 


of fn- We use Cauchy’s residue theorem. We know from 
(47) that fn has a first order pole at wq/c — ie and two 
second order poles at ±ikx, pictured on Fig.[^ From (48) 
and ( [50| ) we see that we have to close the integration path 
(Jordan loop) in the lower half of the complex plane for 
||x|| > 0 and ||x|| +ct > 0 for and /i„ respectively, and 
in the upper half of the plane for ||x||< 0 or||x||+ct <0 
for Qn and respectively. 


It can be checked that the residues of (w) e and (w) e read 


Res(5„(-)e-‘'ll"ll,--ie) ^ 

V C / e—>- 0 + 

Res {gn (•)e-‘'ll’^ll,ifcx) (70^^”^ + 7+^”^ l|x||) , 

Res(g„(-)e-‘'ll’'ll,-ifcx) (70”^”^ + 71”^”^ l|x||) 

Res {hn (•,i)e“'''''''', ^ - ie) G^"\ 

Res (h„(-,t)e-HMI,ifcx) (||x|| +ct)) , 

Res(^/i„(-,t)e“''ll^ll,-ifcx) + 7 r^”^ (l|x|| + ct)j 


(51a) 

(51b) 

(51c) 


(52a) 

(52b) 

(52c) 


where the Gq"^ and 7 )*"^"^ depend on n, as suggested by We give 


the notation. One can see that 


+ (n) -*{n) _ (n) 

7o =7o =7o > 

+ (n) —*(n) (n) 

7l = -7i = 7i • 


in) _{‘ff~\ckx) 


Cto — “ 


(Wq + C^fcl) 

\2-n I 


(53a) 


(n) _ ckx (ifcx) [-i (-3 + n) Wo - (n - 2) c/cx] 
7o — 


4 (iwo + cfcx) 


(„) _ ck'^ {ikxY 


7i = 


4 (iwo + ckx) 


(53b) 

(53c) 


The Fourier transform (50) is thus given for n G {1,2} 
by 
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2m 


= -0(l|x||) 


-fcx||x|| _ ^*in) 


7o 


7i l|x|| +e 




+ 9{\\K\\+ct) (7o^”^-7i^”^(||x||+rf)) + 

-0(-(||x||+ci)) (7(”^+7^”^(||x||+ct))' . 

(54) 

We can then compute the convolution—which we call G„ (||x|| ,t )—with the principal value as prescribed by (27). It 
yields, still for n e {1, 2}, 


G„ (||x|| ,t) = - — 


= g-fcxllx 


Vp- 

*(n) *(n) 


(70"" -71' '' llx 


-Ei (fcx I |x| I) + (kx (I |x| I + ct)) 

7_gfcx||x|| _^g(ic^o-cfcx)tgfcx(||x||+ct)A _c7*(”)te(*‘^“-^'=^)*Ei(fcx(||x|| +ct)) 
fcx V ) 


*(n) 

7i 


+ e'=>^ll"ll [( 7 ^") + 7 ^") ||x||) [-Ei(-fcx ||x||) + e(‘‘"“+"'=^)*Ei(fcx (- ||x|| - ct)) 

(n) 

+ ^ [_e-'=xl|x|| + C 7 ^"^t e(*‘^“+"'=^^*Ei (fcx (-||x|| - ct)) 

kx \ ) 


+ e’ 


'G^ 


-Ei (i^||x||) -bEi (i^ (||x||-bet)) 


(55) 


Now, remember that for n = 3, the integrand in (47) has a (simple) pole at fc = 0. As argued above (48), we can shift 
this singularity away from the real axis 0 —?> —ia to compute the integral (see Fig. [^, before taking the limit a —>■ 0 
at the end. This yields 


/3(l|x||,t) 

2i7r 


= -0(l|x||) 


-fex||x|| _ ^*(3) 


i7o 


7i 


'll) 


-b e 


-b 0 (||x|| -b ct) 
-^(-(l|x|| +ct)) 


giwotg-/cx(||x||+ct) 


( *(3) *(3) 

>7o -7i 


’^IIg(^)]+0(-||x||) [e-'=xllx||(^(3)+^(3)||^||^ 
Ixll Tct)) +e-‘^ll’^llGj,^^' 


jii^oigfcx(||x||+ct) /^^13) 


(3) 


X -b ct 


-[0(||x||)-0(||x||+ct)e-«*] 
UJq 


and 


(56) 


G 3 (||x||,t) [( 


*(3) *(3) 11 

i7o -7i l|x 
*(3) 


-Ei (kx I |x| I) + e(^‘^o-ck.)tj^- ( 11^11 + 


2)_ 7_gfex||x|| _^g(i<^o-cfcx)tgfex(||x||+ct)\ _c7*(3)te('"'«-“'=^)*Ei(A:x(||x|| -bet)) 

kx y / 


-be' 


Ax||> 


( 7 ^^) + 7p^ I |x| 1 ) [-Ei (-kx I |x| I) + e(-«+='=-)‘Ei (fc^ (- ||x|| - ct))' 


v(3) 


-b^ [-e-fcx||x|| _^g(ic^o+c/cx)tgfcx(-||x||-ci)A _^g^(3)^g(i<^o+cfex)tEi(fcx(- ||x|| - ct)) 

kx \ ) 






Wo 


-Ei i— X +Ei i— X +ct 


.Wo 


Wo 


-vp 


t*+oo 


di; 


— e 


ILjQt 


p-\-oo 


drt 


(57) 
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which is an infinite quantity, but this is not a problem as we are interested in (— ||x|| ,t) — C 3 (||x|| ,t), which, as 
we shall see, is finite. Keeping in mind that ||x|| and t are both positive we compute, for n € {1, 2} 


fn (- l|x|| ,<) T fn (||x|| ,t) = 2 i 7 r ||x|| +ct)e^~ 


+9{-\\x\\+ct) 


„-fcx(rf-||x|| 


( 70 ' 

(n 


(76 




Te 


_j_gfex(ct- 

( (n) *(n)\ ( (n) *■ 

7o T7 o ) - ( 7 I 'T7i 

-fex(ct+||x||) _ ^*{ 


■7i 

(n) 

■7i 


|x|| — ct) 


(n) 


( 70 ' ' -7i' 


and, for n = 3, the difference 

Is (- l|x|| ,t) - /a (||x|| ,f) = 2 i 7 r ||x|| + ct)e‘- 

+ 0 (-||x||+ct) e*'""* 


(I |x| I + ct)) + - 7^^ (| |x| | - ct) 

(58) 


^-fcx(ct-||x|| 


/ *(3) 

(^7o 


*(3) 

■7i 


|x|| - ct)) 


_|_gfex(ct-| 


(70 


(3) _ ^(3) 


X — ct 


+e-'=xll 


—e 


Te 


/ (3) ^ *(3) 

(^7o T7o 

-/cx(ct+||x|| 


\ / (3) *(3) 

) - (^7i T7i 


*(3)_ ^,*(3)^1^gfex(ci-llxll) |^^(3)_^(: 


1 


- [1 — (1 — 0 (ct — 

Wo 


^7o - 7i 

1^1 h'l ai“ot 


''"xll — ct)) 

(59) 

The quantities G„ (— ||x|| ,t) =p G„ (||x|| ,t) are not illuminating enough to warrant their explicit writing out here, 
but their expression follows immediately from (55) and (|57[) . It is noteworthy, though, that the contribution to 


G 3 (— I |x| I, t) =F G 3 (I|x|I , t) (near-held) from the last summand on the right-hand side of (57) (the summand featuring 
the two principal value integrals) reads 


1 


X T ct 


-Tgi'^ot jgg . 

Wo ^ Vl-llx|| Tctiy 


(60) 


This particular result is reminiscent of the hndings of Karpov et al. in m on the consequences that restricting 
the spectrum to positive frequencies has on localisation and causality. The cited work m dealt with the simpler 
problem of the free propagation of free massless particles in one-dimensional space. An expression similar to (601 was 


derived. As mentioned elsewhere in this manuscript, we intend to return to the specihc features of the near-held in 
an upcoming work. Also, note that (27) can be rewritten 




|x|| ,i) = ^/n (T ||x|| ,t) - G„ (=F I|x|| ,f), 


(61) 


and the photon wave function is given, according to (19), (46) and (47), by 
25 he 


^_L (x,t) = -i 


34 eometto { 2 'Kf 


r c{x,y) 
8>m2 

l|x|| 


1 

T 

(INK^) - ^ 


(l|x||,t) ^ 

\tt{-) rr(+)l 

(l|x|| ,t) 


f(*) 

9 ‘STT12 

L ^iix|i 



IM 


(INI^^) + ||X||5 

- Ht\ 

(l|x|| ,t) 



(62) 


where it is of course understood that the values of all 
functions are taken at (||x||,t). Remember that 

are the components of ^^2 in the plane perpen- 


(z') 

dicular to x, while Qni is the component of ^^2 in the 
direction of x. 
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Figure 5: See rest of figure and caption below. 


VII. DISCUSSION 


According to ([54| and ([^) —and, in the n = 3 case, 


(56) and (57)—our final resu 


evidently noncausal. 
causality: 


t given by (61) and (62) is 


We identify three sources of non- 


The first one is well-known mil] and is due to the 
fact that, in order to compute the photon wave 
function, we integrated over the physical frequen¬ 
cies of the electromagnetic field, which are posit¬ 
ive. Indeed, the spectrum of the electromagnetic 
field Hamiltonian is the positive real semi-axis and 
is thus bounded from below. In light of Hegerfeldt’s 
theorem, it is then not surprising to observe that 
our result is noncausal. More precisely, the theorem 
teaches that the convolution contributions (55) and 


(57), which would not be featured if the integration 


was carried out over the whole real axis, necessarily 


introduce a noncausality in (61). This feature has 
already been studied in sect. VB[ where the dipole 


approximation was performed, which changes little 
to the discussion and the results. 

• The second one is the presence of the singularity 
at A: = 0 of the integrand (47) for n = 3, that is. 


according to (19), in the near-field zone of emission. 


Terms coming from this singularity are given in (56) 


and (57). As can be seen from (62), this part of the 
emitted field decays as the third power ||x||~^ of 
the inverse distance from the nucleus to the point 
of observation. There is much of interest to say 
about the near-field and this noncausality, and how 
they can be linked to the (Coulomb) longitudinal 
electric field, but we will discuss that in a different 
manuscript. 

• The third and last one is the uncertainty on the po- 
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Figure 5: Square moduli of the contributions to the far- (n = 1) and mid- (n = 2) fields from the poles as given by 
(58) and from the principal value integrals as given by (|55[) as a function of distance from the nucleus, for fixed time 


in the left column, and as a function of time, for fixed distance in the right column. Only the contributions to (58) 
and (55) from the poles at A: = ±ifcx are plotted here (see Fig. for those coming from the pole k = ujo/c), which 


means that Gg"^ is artificially brought to zero. This is for the minimal A • p coupling (sect. VI). Both axes are 


logarithmic on all figures. The boundary of the lightcone is signalled by a vertical line. 


sition of the electron. For instance, at t = 0, when 
the emission starts, the electron is in the 2p level, 
and its wave function is spread on a distance of 
order 2ao. In the light of this uncertainty, it is nat¬ 
ural to expect a spacewise-exponentially decreasing 
tail in the emitted held, with characteristic size of 
order ag. And this is what we indeed obtain if we 
neglect the other sources of noncausality, namely if 


we integrate in (46) over both positive and 
negative frequencies, 

we only focus on the mid- (n = 2) and far- 
(n = 1) held zones of emission, 


in other words, if we consider only expression (54), 


the result is still noncausal. In the previous sen¬ 
tence “noncausal” is understood to mean “not van¬ 
ishing outside the lightcone centred around t = 0 
and the position x = 0 of the Hydrogen nucleus 
(proton)”. Rather, the noncausal terms in ( [54| ) de¬ 
cay exponentially on a distance 3ao/2 (see Fig. [^. 
Hence these terms are nonvanishing outside the 
lightcone centred around t = 0 and x = 0, which 
we understand as being an illustration of the fact 
that the electron emits the photon from its own 
position, which is not fully determined and is only 
exponentially conhned within distances of order og 
around the nucleus, rather than from the position 
X = 0 of the nucleus itself [28]. In |2], Shirokov 
very clearly hints at such exponentially decreasing 
tails outside the lightcone. But among the works of 
his which are available to us (including jUj), none 


presents or even mentions the method he used to 
obtain this feature. Note that the contributions to 
the emitted held from the poles aX k = ±ifcx (which 

I I (n) 

correspond to (58) and (55) with Gg artihcially 
brought to zero) which arise when the dipole ap¬ 
proximation is not performed, are also subject, so 
to say, to Hegerfeldt noncausality. This results in 
the (principal value) noncausal contributions seen 
in Fig. [^ Interestingly, these noncausal terms do 
not decay more strongly with increasing distance 
outside the lightcone than they do inside, as was 
the case for the noncausal terms studied in sect. El 
This can be seen either graphically or by noticing 
that the same asymptotic arguments as developed 
in sect. El do not yield the same result here. 


In an upcoming manuscript we will focus particularly 
on the near-held, where the longitudinal (in the sense of 
Fourier space, see EH) contribution to the electric held 
comes into play. 


The results presented here regarding departure from 
causality are, of course, disturbing but maybe not 
exceedingly disturbing. Indeed, one notices from Figs. 
and that the weight of the wave function present out¬ 
side the light cone is always small. In [30] for instance 
we estimated that for a simpler but related problem, the 
weight of the negative frequencies after a time of the 
order of the lifetime is bounded by 10“^°, which could 
never explain for instance non-locality effects a la Bell 
characterized by violations of l/-\/2 (see footnote 3 in 

EOj). 
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Moreover Sipe showed [5] that, performing the same 
kind of approximations as described in sect. IV (usual 
E • X coupling in the dipole approximation, integration 
over the whole real axis), one can approximate QED by 
a theory in which the photon wave function obeys (com¬ 
plex) Maxwell-type equations and thus admits a causal 
(retarded) Green’s function. This explains why depar¬ 
tures from causality are always small. It is not clear 
at this level whether QED is intrinsically non-causal or 
whether apparent non-causalities result from approxim¬ 
ations performed somewhere in the theoretical develop¬ 
ments. 
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Let us note that, for n € {1,2}, the integrand in (34| 


[24] 


[25] 

[26] 


at Do/c is also artihcial. To compute the integral, though, 
we split the integrand into two meromorphic functions in 
order to be able to use the Jordan lemma. 

M. Abramowitz and 1. Stegun, Handbook of Mathemat¬ 
ical Functions with Formulas, Graphs, and Mathematical 
Tables (Dover, 1964). 

This was also true of the integrals for n = 1, 2 in 

the E ■ X case of the previous sect. 
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Note that it is very easy to show that if, instead of this 
time-dependent perturbation theory, we plug in the ex¬ 
pression from the usual Wigner-Weisskopf approximation 
(which has been shown [20] to provide a very good ap¬ 
proximation to the dynamics of the system (except at 
very short times), up to very small corrections), our res¬ 
ults for the matrix elements of the electric field hold, with 
a grain of salt: loq should be replaced by oiQ-l-aiLS —(i/2) P. 
Therefore, our perturbative treatment works directly 
at short times and indirectly (through the substitution 
which we just specified) at intermediate times. It only 
fails at very long times, where the decay becomes nonex¬ 
ponential |20| . 
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It is only if the dipole approximation is performed that 
the emission can be considered to take place at x = 0. As 
discussed at length in m and m, and also below (|42|, 


actually is an entire function of k, since the singularity 
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the dipole approximation results in divergences coming 
from ultraviolet frequencies. A possibility is then to im¬ 
plement a cutoff on high frequencies (the cutoff being of 
order c/oq) which, of course, will induce a “blurring” of 
the lightcone over distances of order ag, in virtue of the 
Paley-Wiener properties of the Fourier transform. In this 
sense, the noncausality due to the finite size of the decay¬ 
ing electron, first noticed by Shirokov [5] and described 
here in detail, is similar to Hergerfeldt’s noncausality. 
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